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SHALLOW  SNOW  MODEL  FOR  PREDICTING 
VEHICLE  PERFORMANCE 

William  L.  Harrison 


INTRODUCTION 

The  state-of-the-art  of  predicting  vehicle  pcformance 
in  snow  is  still  at  the  pioneering  level.  Although  work 
has  been  done  towards  this  end  for  at  least  25  years, 
there  has  yet  to  appear  in  the  literature  a  prediction 
model  based  on  the  treatment  of  snow  as  a  unique 
material  (Yong  and  Harrison  1978).  It  is  anticipated 
that  this  situation  will  change  in  the  next  few  years 
and  attention  will  be  given  to  the  development  of  a 
proper  snow  failure  model. 

The  purpose  of  this  report  is  to  review  all  past 
efforts  that  have  been  devoted  to  the  prediction  of 
vehicle  performance  in  shallow  snow  layers  and  esta¬ 
blish  the  current  state-of-the-art.  Within  this  context, 
a  performance  model  will  be  recommended  for  use  as 
a  “best  estimate”  of  the  traction  and  resistance  devel¬ 
oped  by  wheeled  and  track-laying  vehicles  in  shallow 
snow. 

In  this  report,  consideration  will  be  given  to  the 
modification  of  traction  and  resistance  forces  by 
vehicle  morphology  and  various  types  of  traction  aids. 
The  snow  cover  considered  will  range  in  strength 
properties  from  those  encountered  in  dry  snow  to 
those  of  slush  layers,  and  the  snow  will  be  presumed 
to  be  deposited  on  a  firm  base.  An  assessment  will 
be  made  of  the  accuracy  of  current  methods  and  a 
model  will  be  based  on  available  test  data  from  past 
exercises. 


HISTORICAL  REVIEW 

The  literature  relating  specifically  to  vehicle 
mobility  through  shallow  snow  (ground  support)  is 


not  extensive.  The  problem  has,  however,  been 
sparingly  addressed  in  the  last  25  years  relative  to 
vehicle  performance  predictions.  It  is  generally  agreed 
upon  in  the  literature  that  the  definition  of  shallow 
snow  depends  on  the  vehicle  of  interest.  The  snow- 
pack  is  said  to  be  shallow  when  the  disturbed  volume 
(pressure  bulb)  beneath  a  traction  element  is  in  con¬ 
tact  with  the  ground. 

Within  this  context  of  vehicle  performance  in 
shallow  snow,  several  different  layer  conditions  must 
be  considered.  The  first  condition,  which  represents 
the  idealized  representation  of  shallow  snow,  is  that 
of  a  homogeneous  layer  of  settled  snow  covering  a 
nondeforming  subsurface.  The  second  condition  is 
that  of  a  packed  snow  layer  over  a  nondeforming 
subsurface. 

Other  surface  conditions  such  as  glare  ice,  slush, 
snow-ice,  and  thawing  soil  are  of  interest  in  this 
study  but  will  be  addressed  later  in  this  report. 

Nuttall  (1957)  described  a  possible  approach  to 
predicting  the  resistance  force  acting  on  vehicles 
traveling  in  shallow  snow  layers.  It  will  become 
apparent  that  many  subsequent  efforts  bear  some 
degree  of  similarity  to  this  proposition.  Nuttall’s 
basic  premise  was  based  on  pressure-density  data 
which  in  general  appearance  is  as  shown  in  Figure  1. 
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Figure  7.  Pressure-density  curves. 


The  relationships  developed  by  Nuttall  were  an 
application  of  the  Bernstein  (1913)  and  Bckkcr  (1952) 
equations  for  the  work  of  compaction: 


R  -SZp  dz 
0 


(D 


Nuttall  expresses  p  in  terms  of  density  and  obtains 


R_=  A  fi_  (1-  V) 

F  «  P  70 


(2) 


where  R  =  snow  resistance  to  motion  (N) 

F  =  vertical  force  (N) 

R/F  =  dimensional  quantity  of  resistance  (N) 
per  applied  load  (N) 
h  =  snow  depth  (m) 

£  =  contact  length  (m) 

pz  =  pressure  at  sinkage  z  (kPa) 

p  =  average  ground  pressure  (kPa) 

T|  =  initial  snow  density  g/cm3 
7c  =  final  snow  density  g/cm3 


Factors  relating  to  the  vehicle 

1 .  Geometry  or/and  configuration, 
i.e.  length,  wheel  base,  dimensional  as¬ 
pects  of  traction  elements,  and  number 
of  propulsion  units  (articulated  configu¬ 
ration,  etc.). 

2.  Weight. 

Factors  relating  to  the  operation  of  the 
vehicles 

3.  Speed,  including  slip. 

4.  j  -  slip  distance  of  a  tractive  ele¬ 
ment  relative  to  the  undisturbed  ma¬ 
terial  during  one  cycle  in  contact  with 
the  material. 

5.  V  =  a  characteristic  speed  related 
to  the  vehicle,  such  as  peripheral  wheel 
speed  relative  to  the  vehicle.  At  a  given 
moment,  all  other  velocities  and  compo¬ 
nents  in  the  system  may  be  expressed  in 
terms  of  ratios  to  this  speed.  One  ratio 
of  particular  importance  in  dynamic 
systems  is  the  slip  ratio  I: 


If  eq  2  were  expressed  in  terms  of  resistance  per 
unit  change  in  volume  and  density  in  terms  of  sink- 
age,  it  would  take  the  following  form: 

R  =  bpzz  (3) 

(since  pt  =F/bi  and  z  =  h[1-(yly0)]  )  where  b  is 
track  or  tire  width  and  h  snow  depth.  Equation  3 
bears  a  close  resemblance  to  the  form  postulated  by 
Harrison  (1975)  which  is  discussed  later. 

Nuttall  and  McGowan  /1 961)  collaborated  on  work 
using  scale  models  of  vehicles  in  snow  to  predict  the 
performance  of  full-sized  vehicles.  The  basic  rela¬ 
tionship  developed  and  used  to  express  performance 
was  DlWvs  the  load  numeric  (W/c  dl),  where  D  is 
drawbar  pull,  d  tire  diameter,  W  vehicle  weight,  and 
c  the  structural  cohesion.  This  relationship  comes 
from  the  dimensional  equations  of  the  dependent 
variables  of  sinkage  z,  vehicle  trim  a,  and  drawbar 
pull  D. 

A  significant  feature  of  the  work  by  Nuttall  and 
McGowan  (1961)  is  that,  when  establishing  dimensional 
equations,  much  thought  was  given  to  the  factors 
relating  to  the  vehicle,  to  the  snow,  and  to  the  system 
as  a  whole.  The  following  listing  from  Nuttall  and 
McGowan  (1961)  is  considered  worthwhile  repeating 
as  most  of  these  factors  must  be  considered  in  this 
study.  Some  items  have  been  updated  or  modified 
by  more  recent  research  efforts. 


/=  VJV 

where  is  the  slip  speed  or  average 
speed  of  the  wheel  contact  area  relative 
to  the  undisturbed  soil  or  snow.  It  is 
to  be  noted  that,  while  the  slip  ratio  is 
usually  conceived  in  terms  of  speed,  it 
also  has  a  simple  geometric  interpretation: 

/=//» 

where  /  is  defined  in  factor  4  and  E  is 
the  distance  from  front  to  rear  of  the 
running  gear  contact  area  along  the  same 
path  over  which  j  is  measured  (usually 
parallel  to  the  average  surface  of  the 
material).  In  terms  of  speeds,  the  slip 
ratio  implies  dynamic  snow  reactions 
arising  from  inertial  and/or  viscous 
effects.  In  terms  of  distances,  the  ratio 
implies  essentially  static  material  reac¬ 
tions,  determined  primarily  by  displace¬ 
ment  but  little  influenced  by  time-rate 
of  displacement.  (While  both  effects 
may  be  implied,  dimensional  analysis 
will  not  tell  you  which.) 

Factors  relating  to  the  snowpack 

6.  h  -  total  depth  of  material  of 
interest.  Boundaries  of  significant 


layers,  and  hence  thicknesses  as  well,  may  be 
specified  by  measurements  from'the  surface. 

7.  c  =  before-collapse,  or  structural  co¬ 
hesion,  of  the  material  at  depth  h.  This 
will,  in  general,  vary  at  a  given  spot  and 
time  from  layer  to  layer,  and  hence  with 
depth. 

8.  =  the  full,  unit  after-collapse,  or 
dynamic  shearing  resistance,  of  the  mate¬ 
rial  originally  at  depth  h.  For  any  depth, 
it  may  be  approximated  by  Coulomb’s 
assumption: 

r0  =  c  +  a  tan  #  (6) 

where  a  is  the  unit  normal  loading  on  the 
shear  plane.  On  this  basis,  may  be  re¬ 
placed  by  c  and  tan  #,  where  these  are, 
respectively,  the  effective  cohesion  and  the 
tangent  of  the  effective  angle  of  internal 
friction  of  the  collapsed  material  origi¬ 
nally  at  depth  h. 

9.  y  =  the  before-collapse  bulk  specific 
weight  of  the  material  at  depth  h. 

10.  p  =  kinematic  viscosity  at  depth  h. 

1 1 .  tan  0  =  slope  of  The  surface  of  the 
material. 

12.  f,  =  coefficient  of  friction  of  soils 
or  snows  at  depth  h  in  the  material. 

13.  B  =  stress-strain  parameter  character¬ 
izing  dynamic  shearing  resistance  in  the  layer. 
Dynamic  shearing  resistance  develops  only 
after  some  consolidation  and  reorientation 
of  the  grains  of  the  collapsed  material.  The 
typical  relationship  in  weak,  loose  materials 
between  unit  shearing  resistance  and  shear 
travel  may  be  approximated  by 

T“*0  (1  -erBt)  (7) 

Because  Bi  must  be  dimensionless,  the  di¬ 
mension  of  the  parameter  B  is  S'1 . 

The  list  might  be  further  lengthened  by 
the  inclusion  of  soil  or  snow  elastic  prop¬ 
erties,  but  the  phenomena  of  interest  regu¬ 
larly  involve  large,  permanent  deformations 
of  the  material,  so  that  elastic  forces  may 
be  neglected  from  the  outset,  at  least  until 
experimental  evidence  of  their  importance 
dictates  otherwise. 

Factors  relating  to  the  system  as  a  whole 

14.  g  =  acceleration  of  gravity. 


Dependent  variables 

15.  i  -  sinkage  of  vehicle. 

16.  a  =  trim  of  vehicle 

17.  D  =  drawbar  pull,  or  measurable 
margin  of  tractive  capacity  over  external 
motion  resistance.  Drawbar  pull  of  a 
given  machine  in  a  given  material  is  largely 
influenced  by  slip  ratio  /  and/or  grouser 
travel  /.  Either  drawbar  pull  or  slippage 
may  be  taken  as  the  independent  variable, 
and  the  remaining  one  as  dependent.  In 
testing,  it  is  convenient  to  control  slippage 
and  to  measure  drawbar  output.  In  prac¬ 
tice,  the  drawbar  load  is  fixed,  at  any 
given  moment,  by  terrain  and  towed  load, 
and  slippage  becomes  the  dependent 
variable. 

Most  of  these  factors  will  appear  in  the  algorithm 
for  computing  gross  tractive  effort  and  snow  resis¬ 
tance  to  motion. 

Blackmon  and  Rula  (1960)  presented  a  compila¬ 
tion  of  tests  conducted  in  Boulder  and  Camp  Hale, 
Colorado;  Fort  Churchill  and  Kapuskasing,  Canada; 
and  Houghton,  Michigan.  Some  of  the  tests  can  be 
considered  as  shallow  snow  tests.  Figure  2  shows 
snow  depth  vs  rut  depth  of  two  Army  vehicles  op¬ 
erating  at  Camp  Hale,  Colorado:  a  M8E2  cargo 
tractor  and  a  M5A4  high-speed  tractor. 

Plots  of  the  performance  of  the  M8E2  and  the 
M5A4  (ground  pressure  p  42.7  and  57.9  kPa,  respec¬ 
tively)  indicate  that  the  compacted  snow  density  be¬ 
neath  the  ruts  remained  the  same  regardless  of  the 
original  snow  depth;  i.e.  rut  depth  =  snow  depth  x 
(7/7’).  Assuming  that  the  virgin  snow  density  7 
was  constant,  the  compacted  snow  density  7 'would 
also  have  to  be  constant  in  order  that  the  slope  of 
the  curves  in  Figure  2  remain  constant.  Figure  3 
(again  from  Blackmon  and  Rula  1960)  shows  the 
relationship  between  maximum  drawbar-pull  and 
snow  depth  for  the  M8E2  and  the  M5A4. 

Figures  2  and  3  indirectly  show  the  relationship 
between  rut  depth  and  resistance  in  shallow  snow 
for  the  M8E2  and  M5A4.  Figure  4  (from  Abele 
and  Parrott  1968)  gives  a  tractive  coefficient  for 
the  same  two  vehicles  of  tan#  =  0.44. 

The  tests  conducted  at  Houghton,  Michigan,  by 
Abele  and  Parrott  (1968)  did  not  include  the  same 
type  of  data  as  did  tests  at  Camp  Hale,  and  no 
assessment  could  be  made  of  the  two  vehicles  most 
likely  to  give  shallow  snow  test  results:  namely  the 
D4  and  D7  tractors  in  approximately  46-cm  snow 
depth.  One  very  interesting  result  which  comes 
from  Table  24  of  Abele  and  Parrott  (1968)  (a 
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Table  1.  Shallow  snow  data-snow  characteristics. 


Location 

r* 

ro 

y  , 

(y/cm3) 

h 

(cm) 

(cm) 

dt 

(cm) 

CJ 

(loules/cm3) 

cu 

(kPa) 

<l> 

(degrees) 

c 

(kPa) 

<> 

(degrees) 

Houghton 

-  7 

0.15 

14 

10 

20.3 

0.30 

0 

18.8 

U 

18.2 

30  |an  75 

-  7 

0.15 

16 

12 

20.3 

0.50 

0 

18.8 

1.5 

18.2 

-  7 

0.15 

14 

11 

20.3 

0.35 

0 

18.8 

1.5 

18.2 

-  7 

0.15 

11 

9 

20.3 

0.29 

0 

18.8 

1.5 

18.2 

-  7 

0.15 

12 

10 

20.3 

0.42 

0 

18.8 

1.5 

18.2 

-  7 

0.15 

14 

10 

25.4 

0.29 

0 

18.8 

1.5 

18.2 

-  7 

0.15 

17 

13 

20.3 

0.36 

0 

18.8 

1.1 

20.2 

-  7 

0.15 

18 

14 

20.3 

0.64 

0 

18.8 

1.1 

20.2 

-  7 

0.15 

19 

IS 

20.3 

0.55 

0 

18.8 

1.1 

20.2 

-  7 

0.15 

18 

14 

25.4 

0.59 

0.27 

17.5 

1.1 

20.2 

-  7 

0.15 

21 

IS 

20.3 

0.49 

0.27 

17.5 

0.89 

20.2 

-  7 

0.15 

21 

16 

20.3 

0.63 

0.27 

17.5 

0.89 

20.2 

-  7 

0.15 

20 

16 

20.3 

0.66 

0.27 

17.5 

0.89 

20.2 

-  7 

0.15 

19 

13 

25.4 

0.42 

0 

20.3 

0.89 

20.2 

-  5 

0.14 

17 

12 

20.3 

0.55 

0 

20.3 

0.4 

22.2 

-  5 

0.14 

16 

12 

20.3 

0.50 

0 

20.3 

0.4 

22.2 

-  5 

0.14 

16 

11 

20.3 

0.36 

0 

20.3 

0.4 

22.2 

-  5 

0.14 

13 

10 

20.3 

0.48 

0 

20.3 

0.4 

22.2 

-  5 

0.14 

14 

11 

20.3 

0.57 

0 

20.3 

0.4 

22.2 

-  5 

0.14 

16 

12 

20.3 

0.58 

0 

20.3 

0.4 

22.2 

-  5 

0.14 

14 

10 

25.4 

0.31 

0 

20.3 

0.45 

22.2 

19  Feb  75 

-12 

0.26 

11 

7 

20.3 

0.60 

0 

22.8 

0.13 

25.9 

-12 

0.26 

11 

8 

20.3 

0.91 

0 

22.8 

0.13 

25.9 

-12 

0.26 

11 

8 

20.3 

0.60 

0 

22.8 

0.13 

25.9 

-12 

0.26 

10 

7 

20.3 

0.70 

0 

22.8 

0.13 

25.9 

-12 

0.26 

11 

6 

25.4 

0.40 

0 

22.8 

0.13 

25.9 

-12 

0.26 

12 

6 

25.4 

0.60 

0 

22.8 

0.13 

25.9 

-  4.5 

0.23 

14 

9 

20.3 

0.28 

0.34 

25.2 

0.69 

24.7 

-  4.5 

0.23 

11 

9 

20.3 

0.58 

0.34 

25.2 

0.69 

24.7 

-  4  5 

0.23 

12 

5 

20.3 

0.38 

0.34 

25.2 

0.69 

24.7 

-AS 

0.23 

9 

6 

25.4 

0.26 

0.34 

25.2 

0.69 

24.7 

3  Mar  75 

-  3 

0.15 

9 

7 

20.3 

0.17 

0.27 

22.8 

-  3 

0.15 

10 

7 

20.3 

0.28 

0.27 

22.8 

-  3 

0.15 

12 

8 

20.3 

0.27 

0.27 

22.8 

-  3 

0.15 

7 

6 

25.4 

0.05 

0.27 

19.4 

-  3 

7 

6 

20.3 

0.14 

0.27 

19.4 

-  3 

8 

6 

20.3 

0.21 

0.27 

19.4 

-  3 

0.30 

12 

7 

20.3 

0.38 

0 

19.8 

-  3 

0.30 

12 

9 

20.3 

0.50 

0 

19.8 

-  3 

0.30 

13 

0 

20.3 

0.32 

0 

19.8 

-  3 

0.30 

16 

10 

25.4 

0.35 

0 

19.8 

4  Feb  76 

11 

8 

15.2 

1.27 

1.0 

21 

10 

7 

15.2 

1.01 

1.0 

20 

10 

8 

15.2 

1.42 

1.0 

18 

8 

7 

15.2 

1.34 

1.0 

19 

5  Feb  76 

9 

7.7 

15.2 

1.61 

1.0 

21 

11 

9.3 

15.2 

1.24 

1.4 

21 

18  Feb  76 

15 

10 

15.2 

0.69 

1.2 

33 

14 

11 

15.2 

0.92 

1.2 

33 

17 

12 

15.2 

0.80 

1.2 

33 

18 

12 

15.2 

0.75 

1.2 

33 

19  Feb  76 

16 

13 

15.2 

1.33 

0.68 

24.7 

16 

13 

15.2 

1.32 

0.68 

24.7 

15 

12 

15.2 

1.39 

0.68 

24.7 

*  Snow  temperature, 
t  Plate  diameter. 


^*^‘*'-*1'  jfc^f'ara jj> 


Table  1  (cont’d).  Shallow  snow  data-snow  characteristics. 


Location 

rc> 

(g/cm3) 

h 

(cm) 

z 

(cm) 

d 

(cm) 

UJ  ca  41 

(joules/cm3)  (kPa)  (degrees) 

c 

(kPa) 

4> 

(degrees) 

Silver  Creek  Bridge  * 

-  5 

0.1 

29.2 

17 

5.0 

14  4+ 

0.44 

20.8 

30  Jan  73 

-  i 

0.1 

29.2 

15 

7.6 

9.00 

0.44 

20.8 

-  5 

0.1 

29.2 

13 

10.2 

5.03 

0.44 

20.8 

-  5 

0.1 

29.2 

9 

12.7 

2.80 

0.44 

20.8 

-  5 

0.1 

29.2 

16 

7.6 

9.10 

0.02 

31.2 

-  5 

0.1 

29.2 

15 

10.2 

5.65 

0.02 

31.2 

-  5 

0.1 

29.2 

12 

12.7 

3.4 

0.02 

31.2 

-  5 

0.12 

25.4 

15 

7.6 

in 

0 

27 

-  5 

0.12 

25.4 

11 

10.2 

1  U.JJ 

4.59 

0 

27 

-  5 

0.12 

25.4 

9 

12.7 

2.30 

0 

27 

-  6 

0.12 

25.4 

19 

5.0 

16.23 

0.15 

28.1 

-  6 

0.1 

25.4 

15 

7.6 

8.39 

0.15 

28.1 

-  6 

0.1 

25.4 

11 

10.2 

3.72 

0.15 

28.1 

-  6 

0.1 

25.4 

11 

12.7 

3.68 

0.15 

28.1 

-  4 

0.16 

23.5 

22 

5.0 

16.35 

0.11 

26.6 

-  4 

0.16 

23.5 

20 

7.6 

7.23 

0.11 

26.6 

-  4 

0.16 

23.5 

21 

10.2 

5.45 

0.11 

26.6 

-  4 

0.16 

23.5 

22 

12.7 

4.03 

0.11 

26.6 

-  6 

0.15 

25.4 

22 

5.0 

13.17 

0 

31.2 

-  6 

0.15 

25.4 

23 

7.6 

7.33 

0 

31.2 

-  6 

0.15 

15.4 

23 

10.2 

4.98 

0 

31.2 

-  6 

0.15 

23 

20 

12.7 

4.17 

0 

31.2 

-  6 

0.15 

23 

21 

12.7 

3.38 

0 

31.2 

-  8 

0.14 

40 

29 

7.6 

12.18 

0.27 

23.5 

Ft.  Greely 

2  Feb  72 

-  7 

0.2 

43 

34 

11 

1.13 

2.4 

32 

3  Feb  72 

-12 

0.2 

51 

48 

11 

1.08 

1.6 

34 

14  Feb  72 

-12 

0.22 

58 

48 

11 

1.05 

2.3 

36 

IS  Feb  72 

-10 

0.22 

51 

39 

11 

0.92 

1.6 

32 

Ft.  Walnwrlght 

23  Feb  72 

-10 

0.2 

61 

58 

11 

1.57 

2.3 

34 

24  Feb  72 

-  8 

0.2 

59 

56 

11 

1.31 

0.75 

38 

28  Feb  72 

-10 

0.21 

74 

63 

11 

1.30 

2.3 

36 

1  Mar  72 

-10 

0.21 

56 

51 

11 

1.62 

1.3 

34 

Camp  Hale,  Colorado 

24  Feb  58  to 

-  8 

0.25 

76  • 

35 

0.92 

7  Mar  58 

-  7 

0.27 

77 

58 

2.26 

-  6 

0.24 

74 

49 

1.99 

-11 

0.27 

71 

52 

2.20 

Wolf  Creek  Rd. 
California 

-1 

0.36 

14 

24 

2.11 

21  March  80 

-1 

0.36 

14 

- 

24 

1.85 

-  1 

0.36 

25.4 

12.9 

24 

1.86 

-  1 

0.38 

22 

9.3 

24 

2.28 

A  Ita,  Utah 

26  March  80 

0 

0.13 

6.3 

24 

2.35 

27  March  80 

-3 

0.13 

24.4 

. 

24 

0.67 

-.3 

0.13 

25.9 

21 

24 

1.08 

30  March  80 

-7 

0.095 

26.7 

19.5 

24 

0.972 

1  April  80 

-1 

0.18 

30.3 

25.5 

24 

1.42 

** 

-  1 

0.18 

22.9 

13.7 

20.32 

0.720 

** 

-1 

0.18 

22.9 

13.7 

20.32 

0.682 

♦  Bennett  (1974) 

}  Large  values  due  to  a  compacted  layer  of  snow  6  to  10  cm  under  new  snow. 
•*  Plate  test 
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0  in  20  30  40  50 

Nominal  Ground  Contact  Prtttur*  (KFo) 

Figure  4.  Compzr  'on  of  nominal  ground  contact  pressure  and 
unit  tract've  effort. 


summary  of  the  Houghton  tests)  is  sup  tort  of  the 
critical  density  theory  of  Gerdel  et  al.  (1  }54),  *  .uuuli  ’ 
and  Finelli  (1955),  and  Harrison  (1975).  The  aft:-' 
traffic  density  measurer  t  of  the  f>  st  pass  gave  a 
value  of  approximately  0.40  for  all  \e<  l  ies  regardless 
of  ground  prt:-.!-e  Th:s  was  true  wi-:, ..  tiy  two  ex¬ 
ceptions:  one  wi'h  «.*e*;«r  values  vUu.  i -actor) 
and  tne  .ie  with  'o.-'er  values  (the  M29  weasel). 
Unfortunately,  th.se  tests  we*e  the  u.'h'  ones  from 
Abelc  and  P<.r:-V,t  (1 9(58)  that  coo  <.!  be  'f  nsidered 
as  shallow  s;.o  .  lest'  . 

A  summars  tit*'  •  csults  of  the:t  (ostc  >'  v,» 
viewpoint  o'  -sniclc  peifo'K.tiice  in  shallow  s  tow  is 
presented  by  S.J.  Knight  as  Appendix  C  in  CRkEL 
Technical  Report  268  (Harrison  1975).  The  conclu¬ 
sions  were  that,  while  the  prediction  techniques  and 
equations  developed  from  these  data  were  suitable 
for  predicting  performance  in  deep  snow,  they  were 
not  applicable  to  the  shallow  .now  problem. 

An  expedient  solution  for  wheeled  vehicles  travel¬ 
ing  through  shallow  snow  was  presented  by  Harrison 
(1 973)  for  predicting  motion  resistance  based  on  snow 
reaching  a  critical  density  of  0.4  Mg  m"3 .  Other  per¬ 
tinent  ■  -.umptions  were  that  the  contact  area  for  a 
wheeh.a  vehicle  would  be  the  projected  area  under 
static  load  for  a  specific  tire  pressure.  These  data 
are  usually  available  from  the  tire  manufacturer. 

The  purpose  of  the  approach  presented  in  Harrison 
(1 973)  was  to  compute  HmiX  with  no  consideration 
given  to  slip  or  interface  velocity. 

To  predict  gross  tractive  effort  H  and  motion  re¬ 
sistance  R,  the  following  equations  were  suggested. 

For  gross  traction: 

H  =  Wp  +  Ac '  (8) 


W  -  vehicle  weight  (or  total  load  on  powered 
wheels,  kg) 

p  =  tan  41  where  0  <  4»  <  0  (angle  of  inter¬ 
face  friction) 

A  -  nominal  contact  area  (static  tire  print  for 
wheels,  mJ) 

c'-c  — . _  (adhesion)  0  <  c'<  c  (kPaJ) 
tam> 

For  sinkage: 

*  =  (i-7rr)*  0) 


h  =  snow  depth  (m) 

7  =  equivalent  water  content  of  snow 
(density  in  Mg  m‘3). 

For  compaction  resistance: 


Rc  =  0.38  6p°-625  z1,975  atn  =  1.6 
b  =  track  or  tire  width 
k  -  plzn 

p  =  contact  pressure 
For  bulldozing  resistance: 

Rb  =  (M«  +  yz7Ny*+czN*J2b 

R  =  yz2Ny  +  czNti  (N/m) 

7  =  snow  density  (Mg  m-3 ) 
z  =  sinkage  (m) 


where 


wW®' 


c  =  snow  cohesion  (kPa) 

"N"  factors  =  f(0;p)  (12) 

where  0  =  snow  internal  friction  angle 

p  =  (90°  +  approach  angle  of  track)  or  (90° 

+  arc  cos  for  wheeled  vehicles] 

d  =  tire  diameter  (m) 

Equation  10  for  compaction  resistance  is  a  modifi¬ 
cation  of  the  deep  snow  tracked  vehicle  equation  used 
by  the  Land  Locomotion  Laboratory  (Harrison  and 
Czako  1961).  It  was  based  on  a  selected  value  of/7  =  = 

1.6  which  was  representative  of  a  number  of  experi¬ 
ments  in  deep  snow. 

Equation  1 1  for  bulldozing  resistance  was  adapted 
from  the  work  of  Reece  (1965)  and  Hettiaratchi  (1968). 
The  curves  for  determining  the  "N"  factors  are  con¬ 
tained  in  Hettiaratchi’s  work  but  are  not  considered 
relevant  for  inclusion  or  further  discussion. 

Liston  (1974a,  1974b)  presented  an  interesting 
approach  that  uses  the  sinkage  of  a  plate  which  loads 
the  snow  cover  to  the  same  pressure  as  would  a  tire 
having  a  specified  diameter,  width  and  inflation  pres¬ 
sure.  An  abridged  version  of  his  approach  follows. 

It  is  assumed  that  the  maximum  shear  strength  of 
snow  can  be  described  by  the  Mohr-Coulomb  equation: 

r  =  c  +  p  tan  0  (13) 

where  r  =  shear  stress  (psi) 
r  =  shear  stress  (psi) 
c  =  apparent  cohesion  (psi) 
p  =  normal  pressure  (psi) 

0  =  angle  of  internal  friction. 

In  determining  the  magnitude  of  c  and  0,  a  device 
having  a  rubber  surface  (rather  than  a  roughened 
surface)  is  used  to  impose  a  failure  surface  within 
the  snow  mass.  The  more  sophisticated  equation  of 
Hanamoto  and  Janosi  (1961),  which  states  that 

t  =  <c  +  p  tan  0)  (1  -  e~^k),  (14) 

was  not  used  because  it  was  found  that  the  magni¬ 
tude  of  k  was  small,  on  the  order  of  0.02  m,  so  that 
maximum  shear  strength  was  developed  along  most 
of  the  tire/snow  interface. 

The  contact  area  was  assumed  to  consist  of  two 
sections:  1)  an  area,4p  beginning  at  the  uppermost 
contact  point  between  tne  wheel  and  snow  and  ending 
at  the  bottom  of  the  rut  (see  Figs.  5  and  6),  and 
2)  an  area/4  consisting  of  the  flat  contact  patch  at 


the  bottom  of  the  tire.  Ap  is  identified  as  the  peri¬ 
pheral  area  and  A  as  the  contact  area. 

In  order  to  compute  the  peripheral  area  Ap,  it  is 
noted  that 


dAp  =  br  d 0 


as  shown  in  Figure  5,  and 
/«0 


J  o 


dO  =  tr  0Q. 


0O  =  cos"1  (  r  •) 


A  =  br  cos-1  (  7  )  (18) 

In  computing  the  contact  area  (A),  it  is  assumed 
that 

W~P\A  (19) 


- 

-*Vd  9 

Figure  5.  Determination  of  wheel / 
snow  interface  area  (from  Liston 
1974a,  1974b). 


t\4  vj  ?  n< 


Figure  6.  Volume  change  produced  by  passage  of 
wheel  and  identification  of  contact  areas. 


where 


in  which 


p,  =  inflation  pressure  (psi) 

A  =  contact  area  (in.2) 

Equation  19  is  inaccurate  as  it  does  not  account  for 
the  support  provided  by  the  tire  carcass.  The  total 
contact  area  A  t  per  wheel  is  taken  as 

A=(A-^A)  (20) 

and  gross  traction  is  found  from 

H  =  (/1p  +  A)  c  +  W  tan0.  (21) 

In  determining  motion  resistance,  it  was  assumed 
that  the  work  expended  in  compacting  the  snow  was 
the  single  source  of  resistance  offered  by  the  snow. 
The  following  additional  assumptions  were  made: 

1.  As  shown  in  Figure  6,  the  only  dimensional 
change  that  occurs  when  the  snow  is  compacted  from 
its  original  volume  U,  to  its  final  volume  V 2  is  in  the 
snow  depth  h. 

2.  The  pressure  along  the  interface  of  the  wheel 
and  snow  is  zero  at  the  point  of  contact  and  is  equal 
to  the  inflation  pressure  at  the  lowest  point  of  con¬ 
tact. 

3.  The  variation  in  pressure  along  the  whecl/snow 
interface  is  hyperbolic  in  form.  The  pressure  distribu¬ 
tion  is  depicted  in  Figure  7. 

The  work  expended  in  compacting  the  snow  from 
k',  to  kj  is 


The  assumption  that  the  pressure  p  varies  hyperboli- 
cally  leads  to  the  following  relationship: 

K 

P=P\y  P3) 


Figure  7.  Assumed  pressure  distribu¬ 
tion  along  wheel /snow  interface. 


Pj  =  the  inflation  pressure 
V j  =  bh]i  =  the  initial  volume  of  a  snow  mass 
k'j  =  bh2i  =  the  final  volume  of  the  same  snow 
mass 

Introducing  24  into  eq  23  and  completing  the 
integration  gives  the  following  expression  for  com¬ 
paction  energy  E : 

E  =  P|  V.  In  — , 

12  l/, 

If  /?r  is  defined  as  the  force  resisting  motion  which 
is  related  to  compaction,  then  work  can  be  expressed 
by  multiplying  Rc  by  2.  Also,  the  volume  can  be 
expressed  as  bhxt  Substituting  the  above  expressions 
for  E  and  V  gives  the  equation 
h2 

^c=P^2inf^'  W 

The  total  motion  resistance  Rx  is  obtained  by  adding 
the  compaction  resistance  to  the  hard  surface  rolling 
resistance  RH .  In  order  to  predict  net  tractive  effort 
Hn,  which  is  the  parameter  of  interest,  the  motion 
resistance  is  subtracted  from  the  gross  traction: 

-2  I*.,  +  *„)  <26> 

in  which  /?cl  and  Rc2  are  the  resistances  of  the 
front  and  rear  wheels,  respectively.  Values  of  h2 
are  obtained  from  plate  tests  where  the  contact 
pressure  of  the  plate  and  the  tire  are  approximately 
equal. 

Harrison  (1975)  presented  a  method  for  deter¬ 
mining  the  motion  resistance  of  vehicles  in  shallow 
snow  based  on  the  determination  of  a  unit  energy 
constant.  The  gross  tractive  effort  is  found  by  the 
use  of  a  now-familiar  expression: 

H  =  Aca+W  tan  <p  (27) 

where  A  is  the  total  contact  area  of  the  traction 
elements  and  W  is  the  vehicle  weight. 

Motion  resistance  over  shallow  snow  is  defined  as 
the  work  performed  by  the  vehicle  in  compacting 
undisturbed  snow  when  forming  wheel  ruts.  Figure 
8  shows  the  stress-deformation  curve  of  a  column  of 
snow  under  simple  compression. 


The  equation  for  work  £  is  stated  simply  as 
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Figure  8.  Column  of  mow  under 
simple  compression. 

which  is  basically  the  same  as  those  of  Liston  (1974), 
Bekker  (1956),  Nuttall  (1957),  etc.  The  characteris¬ 
tics  of  the  curve  differ  from  simple  compression  in  a 
plastic  medium  in  that  there  is  compaction  throughout 
deformation  z  rather  than  plastic  flow  as  described 
by  the  typical  soil  model.  In  addition,  the  deforma¬ 
tion  z  approaches  a  maximum  valuez^as  maximum 
compaction  (critical  density)  is  approached  and  the 
o  -  z  curve  becomes  asymptotic.  Critical  .density  is 
assumed  to  be  approximately  0.5  Mgm"3.  This 
value  of  y  is  considered  acceptable  for  two  reasons: 
first,  the  additional  deformation  between  y  =  0.5 
and  7  =  1.0  is  minute  and  the  increase  in  o  required 
is  quite  large,  and  secondly,  in  the  unconfined  slate 
that  exists,  plastic  flow  will  occur  at  this  density 
rather  than  a  decrease  in  volume  (Kinoshita  ar  d 
Akitaya  1970). 

The  curves  in  Figure  9  illustrate  the  use  of 
The  effective  sinkage  (ze)  is  dependent  on  the  depth 
at  which  Vi  7max  occurs  and  is  easily  computed 
from  the  following  equation  for  predicting  sinkage 
(Harrison  1973): 


(29) 


where  y  '  is  the  selected  value  of  post-compaction 
density.  In  practice,  the  measured  curves  are  final- 


Figure  9.  Plate  sinkage  tests 
(Harrison  1975). 


zf  line  and 


ized  (when  necessary)  by  locating  th 
constructing  the  finished  curve  so  that  it  becomes 
asymptotic  with  ze  (Fig.  9). 

The  area  beneath  the  load-sinkage  curve  in  Figure  9 
represents  the  total  work  E  expended  in  compacting 
the  snow  from  z  3  0  to  z  =  2e>  To  obtain  the  work 
of  compaction  per  unit  volume  o>  the  total  work  £ 
is  divided  by  the  snow  depth  h  and  the  area  of  the  cir¬ 
cular  plate  A  p  used  to  obtain  the  curves: 


~  E 


(30) 


HA. 


To  minimize  the  effect  of  size,  a  minimum  plate 
diameter  of  20  cm  is  used. 

The  resistance  to  motion  from  snow  compaction 
is  given  by  the  equation 


Rc  =  2b<wh 


(31) 


where  b  is  the  effective  wheel  width  (or  track  width). 

Bekker  (1976)  presented  a  method  of  determining 
rolling  resistance  based  upon  the  shallow  snow  tests 
conducted  by  Harrison  (1975).  The  method  of  so¬ 
lution  combined  the  earlier  efforts  of  Bekker  and 
Semonin  (1975)  for  determining  the  rolling  resistance 
due  to  tire  flexing  and  those  of  Bekker  (1973)  for 
determining  the  resistance  to  motion  by  vehicles 
traveling  in  a  layered  medium. 

To  determine  the  snow  resistance  to  motion, 
Bekker  utilizes  a  multi-layered  approach  based  on 
the  characteristics  of  plate  sinkage  curves  in  shallow 
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s. 


snow.  As  shown  in  Figure  10,  the  p  -z  curve  is 
plotted  in  log-log  form  and  the  “zones"  are  selected 
dependent  on  obvious  changes  in  slope  “n."  Each 
zone  produces  a  set  of  Ac,  k^,  and  n  values  by  which 
the  resistance  for  that  zone  can  be  calculated. 

ror  "zone  1  ”  where  n  is  generally  less  than  3.0 
or  for  any  zone  where  n  <  3.0,  the  resistance  is 
determined  by 

*C =  MIlI!  (32) 

n  +  1 

where  z  is  determined  from  the  expression: 


From  Bekker  (1976) 

[—■Li 

=  0.3  bQk  (3  -n)  s/D  z,  *  (34) 

Whenever  there  is  a  vdseren  >  3.0,  the  “full 
solution  of  the  integral  W  -  b0  f*  px  dx  must  be 
used"  (Bekker  1976)  so  that  Z° 


a  =  (1  -  0.509  n.  +  0.222  nf  -  0.052  nf  + 

+  0,005  nf  )  (36) 


and 

b.  =  (0.25  -  0.26  n.  +  0.1 37  nf  -  0.028  nf  )(37) 

Rc  is  then  the  sum  of  /?(,  /?(|...  /?j  and  the  total 
resistance  for  a  wheeled  vehicle  is 

R  =2(Rc  +  Rt)  +  ny/Rt  (38) 

where  r?w  is  the  number  of  trailing  wheels,  i.c.  2  for 
a  4  x  4,  6  for  an  8  x  8,  etc. 

Nuttall  et  al.  (1975)  developed  an  ad-hoc  shallow 
snow  model  for  wheels  and  tracks.  The  model  is 
used  as  part  of  the  Army  Mobility  Model  to  predict 
traction  and  resistance  of  a  snow  cover  deposited  on 
frozen  soil. 

The  model  computes  maximum  traction  H  and 
resistance  R.  Net  traction  is  determined  as  the 
difference  between  H  and  R: 

H  =  cA  +  Wn  (39) 

where 

c  =  interface  shearing. resistance  (psi) 

A  -  contact  area  (in.  2) 

W  =  vehicle  weight  (lb) 

H  =  tangent  of  the  angle  of  interface  shearing 
resistance. 

For  wheels  the  following  expression  is  used: 

=12.  "JL  iL  (4°) 

W  nt  d  8 

where 

?  ~  ^  ^6 d  -  52  (in.) 
d  -  tire  diameter  (undeflected)  (in.) 


Figure  1 0.  Pressure-slnkage  graph  of  load-slnkage  function  In  shallow  snow. 
Snow  depth  Is  approximately  9  cm,  temperature  Is  -  1°C,  and  plate  diameter 
about  20.3  cm  (adapted  from  Harrison 's  1975  data). 
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5  =  tire  deflection  (in.) 

/?t  =  maximum  number  of  tires  per  axle 


«a  =  number  of  axles 


b  =  tire  section  width  (in.] 
h  =  snow  depth  (in.) 

7  =  snow  specific  weight 


For  tracks,  Nuttall  (1957)  uses  a  simpler  expression 


JL=  (L  -  0.15) 

W  \L  / 


(41) 


where  L  =  overall  track  length  (in.) 

The  relationships  were  developed  by  Nuttall  from 
a  compilation  of  experimental  test  data  available  in 
the  literature. 

Table  2  gives  a  comparison  of  predictions  by  four 
models.  These  predictions  were  made  from  past  data 
relative  to  vehicle  tests  in  shallow  snow  as  presented 
in  Liston  (1974a,  1974b),  Harrison  (1975)  and  Nuttall 
(1975).  The  prediction  by  the  Bekker  model  was 
furnished  the  author  by  private  communication  with 
Dr.  Bekker. 

The  high  value  produced  by  the  Liston  model  for 
the  Hanover  tests  results  from  the  model’s  character¬ 
istic  of  inflation  pressure  dependency. 

It  can  easily  be  shown  that  the  Liston  (1974), 
Bekker  (1976),  and  Nuttall  (1957)  models  (since  all 
are  compaction  energy  based)  can  be  resolved  into 
eq  31 ,  /?c  =  2bwh,  which  is  the  Harrison  (1 975)  model 
with  the  parameter  u,  the  work  of  compaction  per 
unit  volume,  as  the  basic  transfer  mechanism. 


selected  on  review  of  the  past  work,  which  must  be 
considered  as  the  state-of-the-art,  then  a  number  of 
guidelines  or  goals  should  be  contemplated. 

Before  facing  this  selection  process,  let  us  review 
some  common  aspects  of  the  past  work.  All  methods 
discussed  in  the  review  proposed  the  use  of  some  form 
of  the  Mo'  -Coulomb  relationship  for  predicting  the 
tractive  force  developed  by  a  vehicle.  Further,  most 
methods  proposed  that  the  energy  dissipated  in  making 
the  rut  in  snow  be  considered  as  the  resistance  to 
motion  caused  by  compacting  the  snow.  The  most 
apparent  weakness  of  all  the  methods  is  lack  of  vali¬ 
dation  by  field  tests. 

The  guidelines  for  model  selection  under  the  above- 
stated  circumstances  should  therefore  be  based  on  the 
following: 

1.  The  equations  for  traction  and  snow  resistance 
to  motion  should  clearly  address  the  problem.  These 
equations  should  be  sufficiently  fundamental  that 
future  modification  is  possible. 

2.  The  parameters  which  reflect  snow  strength 
characteristics  should  facilitate  prediction  from  basic 
snowpack  properties  and  be  obtainable  by  field 
measurement  without  great  difficulty. 

3.  Because  of  the  lack  of  a  data  bank  on  shallow 
snow-vehicle  performance  parameters,  any  proposed 
model  must  be  considered  provisional  at  best. 


Traction 

The  prediction  of  gross  tractive  effort  is  based  on 
the  Mohr-Coulomb  expression: 


MODEL  SELECTION 

If  the  model  to  be  proposed  in  this  report  is  to  be 


r  =  a  tan$  +  c 

where  r  is  shearing  stress  (kPa),  a  is  normal  stress  (kPa), 


Table  2.  Comparison  of  prediction  models. 


Comauted  values  of  R  INI* 


Measured  values  INI 


Vehicle 

Location 

Liston 

Harrison 

1975 

Nuttall 

1975 

Becker 

K 

rh 

M 1 5 1 A 1 

Silver  Creek  Bridge 
(30  Jan  73) 

2594 

2280 

1584 

2520t 

373t 

2146 

Chevrolet 

Carryall 

Hanover,  N.  H. 

(18  Feb  75) 

1628 

636 

609 

540 

1 2374 

657t 

580 

M151A1 

Houghton,  Mich. 

(18  Feb  76) 

898 

884 

1365 

1490** 

636** 

853 

M151A1 

Houghton,  Mich. 

(30  )an  75) 

862 

643 

676 

1365** 

636** 

729 

*  Compaction  resistance  Rfi  =  R  -  Rj 
t  Measured  by  deceleration  method 
**  Measured  by  towing 
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and  c  and  <p  are  internal  shearing  resistance  and  the 
associated  angle  of  shearing  resistance,  respectively. 

As  used  in  most  vehicle  performance  models,  this 
equation  is  written  in  the  form  for  predicting  gross 
tractive  force  H : 

H  =  Ac+Wxm4>  (43) 

for  vehicles  with  chains  or  grousers.  For  vehicles  with 
pneumatic  tires  or  band  tracks  with  no  grousers, 

H  -  A  ca  +  W  tan  </>  (44) 

where  A  and  W  are  the  contact  area  and  vehicle  weight, 
and  ca  and  0  are  the  interface  shearing  resistance  and 
associated  angle  of  interface  shearing  resistance,  re¬ 
spectively.  The  prediction  of  gross  traction  H  will  be 
in  terms  of  maximum  sustained  force  rather  than  a 
traction-slip  function. 

To  accommodate  traction  aids  in  the  model,  increased 
traction  will  be  reflected  by  multiplying  the  results  of 
a  standard  tire  by  coefficients  obtained  from  experi¬ 
mental  tests  (Table  3)  of  Liston  (1977). 

The  determination  of  total  traction  developed  will 
be  a  summation  of  the  traction  developed  by  each  trac¬ 
tion  element. 

Resistance 

The  energetics  relationship  (Harrison  1975)  which 
states 

=  2  buh  (31 ) 

is  proposed  as  the  resistance  model.  The  relationship 
between  w,  snow  temperature  T,  and  snow  density  7 
shown  in  Figure  1 1  requires  further  validating  by  a 
more  complete  data  bank  than  shown  in  Tables  1  and 
4,  on  which  Figure ,1 1  was  based. 

Slush  and  thawing  soils 

An  ad-hoc  model  is  suggested  for  slush  and  thawing 
soils  until  a  data  bank  has  been  established.  These  two 


materials,  due  to  their  characteristics  of  flowing  rather 
than  compacting  when  disturbed,  will  be  treated  as  sur¬ 
charge  models  having  hydrostatic  properties  relative 
to  the  forces  required  to  cause  flow.  It  is  assumed  that 
the  traction  is  obtained  from  the  firm  supporting  layer. 
The  base  properties  for  slush  will  be  packed  snow  at 
-1°C  (ice),  and  a  layer  having  a  high  degree  of  unfrozen 
water  for  thawing  soil.  The  resistance  force  is  deter¬ 
mined  as  follows: 

Rc  =  JlL  (45) 

2 

where  7  is  the  density  of  the  layer  and  h  is  the  layer 
thickness.  Values  of  7  are  as  follows: 

slushi-7  =  0.75  g/cm3  (750  kg/m3) 
thawing  soil  layers:  7  =  1760  to  2100  kg/m3 

with  loams  on  the  lower  end  and  sands  at  the  higher 
values.  This  relationship  is  loosely  based  on  the  un¬ 
likely  premise  that  there  will  be  no  sinkage  or  break¬ 
through  in  the  supporting  layer. 

Ice,  hard-packed  snow,  packed  snow  (new) 

These  surface  conditions  are  assumed  to  have  neg¬ 
ligible  motion  resistance.  The  emphasis  in  predicting 
performance  is  placed  on  the  following  equations: 

H-  W  tamp  +  Ac  (43) 

or 

H  =  W  tan4>  +/4ca  (44) 

where  p  (tan  <p  or  tan4>)  in  this  context  is  determined 
from  experimental  data.  Figures  12  and  13  show  values 
of  p  tabulated  from  a  large  number  of  experiments, 
as  functions  of  bearing  capacity  and  snow  surface  temp¬ 
erature.  Figure  14  gives  values  of  p  for  ice  relative  to 
surface  temperature.  A  summary  of  Figures  12, 13, 
and  14  is  given  in  Table  5. 


Table  3.  Traction  aid  data 


Pressure  Drawbar  pull  Percent  increase  over  Traction  aid 


Tire 

(kPa) 

(N)_.  . 

NDCC  standard 

coefficient 

NDCC  standard 

55 

2420 

0 

1 

NDCC  standard 

103 

3372 

60 

1.60 

with  chains 

NDCC  radial 

55 

2535 

25 

1.25 

Bias-ply  snow 

55 

3056 

35 

1.35 

Bias-ply  snow 

55 

3252 

45 

1.45 

with  studs 
Summer  radial 
Sand  tire 


55 

55 


3016 

2304 


45 

25 


1.45 

1.25 


c?  -s  -/a  -a  -so 
mmwnt/mr  rev 

Figure  1  1.  w  versus  temperature  and  density. 


100  200  300 

Bparing  Capacity  (lb) 

Figure  12.  p  versus  bearing  capacity-disturbed  snow. 


10  20  30 

Tamparotura  (*C) 

Figure  13.  p  versus  temperature-disturbed  snow. 


zmA-jtsmqgfK 


.  ....  " 


Table  4.  Shallow  snow  data  -  vehicle  characteristics. 


Vehicle 

h 

d 

b 

W 

P 

D 

R 

Location 

type 

(cm) 

(cm) 

(cm) 

(Kg) 

kPa  (psl) 

(N) 

(N) 

Houghton 

30  Jan  75 

M151 

IS 

75 

20 

1,410 

137  (20) 

1,560 

1,312 

13 

75 

21 

1,410 

103  (15) 

1,640 

1,267 

IS 

75 

20 

1,410 

103  (15) 

1,470 

1,036 

13 

75 

22 

1,410 

55  (  8) 

2,090 

1,365 

14 

75 

21 

1,410 

55  (  8) 

1,870 

1,352 

19  Feb  75 

M151 

8 

75 

22 

1,4t0 

103  (15) 

1,380 

720 

16 

75 

1,410 

103  (15) 

1,468 

8 

75 

23 

1,410 

55  (  8) 

1,650 

854 

3  Mar  75 

M151 

12 

75 

21 

1,410 

103  (15) 

2,490 

1,316 

12.7 

75 

22 

1,410 

55  (  8) 

2,220 

1,556 

12  Feb  75 

M151 

14 

75 

22 

1,410 

103  (15) 

2,000 

1,321 

15 

75 

24 

1,410 

SS  (  8) 

2,220 

1,374 

4  Feb  76 

ENGESA 

10 

7751 

172  (25) 

3,430 

5  Feb  76 

215-ton 

14 

775 1 

172  (25) 

4,400 

18  Feb  76 

M151A1 

18 

353* 

1,560 

55  (  8) 

1,490 

14 

353 1 

1,560 

55  |  8) 

1,490 

19  Feb  76 

ENGESA 

215-ton 

16 

775* 

172  (25) 

3,700 

Ft.  Greely 

2  Feb  72 

M113A1 

43 

280* 

38 

9,500 

43.8  (6.3) 

34,000 

3  Feb  72 

SI 

280* 

38 

9,500 

43.8  (6.3) 

39,000 

14  Feb  72 

58 

280* 

38 

9,500 

43.8  (6.3) 

38,000 

15  Feb  72 

51 

280* 

38 

9,500 

43.8  (6.3) 

40,000 

23  Feb  72 

61 

280* 

38 

9,500 

43.8  (6.3) 

32,000 

24  Feb  72 

58 

280* 

38 

9,500 

43.8(6.3) 

34,000 

28  Feb  72 

74 

280* 

38 

9,500 

43.8  (6.3) 

29,000 

1  Mar  72 

56 

280* 

38 

9,500 

43.8  (6.3) 

31,000 

S/A*r  Cree*  Brldgt 

30  Jan  73 

M151A1 

24 

75 

19 

1,452 

55  (  8) 

2,320 

1,214 

24 

75 

19 

1,452 

103  (20) 

1,630 

2,520 

27  Jan  73 

M151A1 

24 

75 

19 

1,452 

103  (20) 

2,380 

1,970 

31  Jan  73 

M34 

1T7 

29 

7,720 

241  (35) 

3,400 

790 

Camp  Hale,  Colorado 

24  Feb  58  to 

M29 

76 

198* 

51 

1,950 

9.6  (1.4) 

5,100 

5,700s 

7  Mar  58 

M5A4 

77 

300* 

31 

9,350 

58  (8.4) 

22,000 

31,808 

M59 

74 

235’ 

53 

17,600 

48  (7.0) 

56,000 

29,100 

M8EZ 

71 

406* 

53 

18,900 

41  (6.0) 

58,000 

35,500 

1  Actual  tire  contact  area  (cm1 ) 

*  Track  length  (cm). 

*  Vehicle  dragging. 


Table  5.  Mechanical  properties  for  ice  and  packed  snow. 


Parameter 

tcc 

Packed  snow 

tan 

0.07 

0.23 

tan* 

0.03 

0.14 

e  (kPa) 

0 

0.3 

ea  (kPa) 

0 

0.6 

15 


T»mp«rotur#  (*F) 


Figure  14.  p  versus  temperature  -  ice. 


Within  the  temperature  ranges  of  the  experiments, 
the  data  scatter  was  significant  with  no  definite  p  vs 
temperature  relationships  apparent.  If  data  were  avail¬ 
able  for  temperature  ranges  well  below  0°F,  perhaps 
some  trends  would  have  been  evident. 

The  contact  area  A  for  pneumatic  tires  will  be  de¬ 
termined  by  the  tire  inflation  pressure  as  follows: 

A  =  JL  (45) 

where  W  is  the  wheel  load  and  P  is  the  tire  inflation 
pressure,  and  by  the  width-length  product  for  tracked 
vehicles: 

A  =  268  (46) 

where  6  and  8  are  the  nominal  width  and  length  meas¬ 
urements. 

Recommendation  for  use  of  eq  45  is  based  on  Table  6 
which  shows  values  of  measured  tire  contact  areas 
divided  by  wheel  load  vs  predictions  using  this  equation. 
The  average  error  is  1 0%  for  eq  45  over  a  number  of 
inflation  pressures  and  tire  structural  characteristics. 

A  10%  average  error  is  quite  acceptable  at  this  time. 

River  and  lake  ice 

The  safe  crossing  thickness  of  river  or  lake  ice,  i.e. 
freshwater  ice,  can  be  determined  from  the  following 
equation  (Johnson  1979): 


where  h  =  ice  thickness,  Cr  =  vehicle  characteristic 
constant  and  ot  =  tensile  strength.  A  reasonable  ten¬ 
sile  strength  for  ice  during  winter  conditions  is  980 
kPa  (Nevcl  1978),  and  in  springtime  this  value  deteri¬ 


orates  to  100  kPa. 

Johnson  (1979)  gives  the  following  values  of  Cc  for 
the  Ml  51 ,  M35A2,  Ml  1 3  and  M60A1  Army  vehicles. 
All  values  are  for  highway  or  combat  loads:  Ml  51  -8.1 
kPa,  M35A2-22.4  kPa,  Ml  13-23.0  kPa,  M60A1-51.7 
kPa.  Thus,  for  example,  during  winter  conditions,  the 
safe  thickness  of  ice  for  the  vehicles  as  listed  will  be: 

Ml  55-0.1 3  m,  M35A2-0.36  m,  M 1 1 3-0.37m, 

M60A1  -0.83  m. 

Model  use 

The  shallow  snow  model  consists  of  the  relationships 
given  eq  42-47.  Where  possible,  values  of  surface 
strength  parameters  have  been,  or  will  be,  furnished 
from  experimental  tests.  When  this  is  not  possible, 
the  values  will  be  predicted.  Toward  this  end,  a  pre¬ 
liminary  study  was  made  (Berger  in  prep.).  These 
studies  will  be  later  expanded  in  scope  to  include  the 
forecasting  of  all  parameters  including  packed  snow, 
slush,  thawing  soils  and  ice. 

Values  of  H  and  R  are  required  by  the  Army  Mobil¬ 
ity  Model  (Nuttall  et  al.  1975).  The  decision-making 
process  for  obtaining  these  values  by  using  eq  42-46 
(and  the  associated  tables)  is  in  fact  the  shallow  surface 
layer  portion  of  the  eventual  "Cold  Regions  Mobility 
Model.”  This  is  intended  for  use  as  a  submodel  to  the 
Army  Mobility  Model  and  associated  models.  The 
following  flow  diagrams  and  examples  present  an  al¬ 
gorithm  for  the  decision-making  processes  for  deter¬ 
mining//  and./?. 

Factors  to  be  considered  in  the  algorithm  are: 

1.  Number  and  arrangement  of  traction  elements 
(i.e.  driven,  towed,  duals). 

2.  Characteristics  of  traction  elements: 

(i.e.  smooth  "rubber,"  grousered,  inflation 
pressure,  chains,  studs,  ply  design)  (see  Tablc3). 

3.  Surface  material:  (sec  Tables  5  and  7). 
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Table  6.  Tire  contact  area  (in.^). 


Measured 

Theoretical 

A 

X  Error 

52.6 

S3.9 

1.3 

2.5 

45.2 

49.3 

4.1 

9.1 

46.1 

49.3 

3.2 

6.9 

50.3 

53.9 

3.6 

7.2 

49.5 

49.3 

+0.2 

0.4 

52.6 

53.9 

1.3 

2.5 

45.0 

49.3 

4.3 

9.6 

49.5 

53.9 

4.4 

8.9 

44.0 

49.3 

5.3 

12,0 

48.5 

53.9 

5.4 

11.1 

50.3 

53.9 

3.6 

7.2 

50.9 

S3.9 

3.0 

5.9 

54.2 

53.9 

+0.3 

0.6 

49.2 

53.9 

4.7 

9.6 

49.3 

53.9 

4.6 

9.3 

46.2 

53.9 

7.7 

16.7 

41.5 

49.3 

7.8 

18.8 

48.1 

53.9 

5.8 

12.1 

47.0 

53.9 

6.9 

14.7 

45.7 

53.9 

8.2 

17.9 

46.8 

53.9 

7.1 

15.2 

45.1 

53.9 

8.8 

19.5 

48.0 

53.9 

5.9 

12.3 

Avg. 

10.0% 

Actual  programming  of  the  model  has  not  yet  been 
undertaken.  A  description  of  the  model  can  be  divided 
into  two  portions,  the  input  data  structure  and  the 
various  subprograms  which  calculate  traction  and 
motion  resistance. 

An  explanation  of  the  data  structure  involves  a 
description  of  all  the  pertinent  variables  and  how  they 
arc  input  into  the  model.  The  computer  model  breaks 
the  vehicle  down  into  a  series  of  units.  Each  transverse 
set  of  axles  or  set  of  tracks  is  considered  a  unit.  The 
variables  pertaining  to  the  geometry  of  wheeled  and 
tracked  vehicles  are  shown  in  Figures  1 5  and  1 6  re¬ 
spectively.  Data  for  particular  tire  or  track  character¬ 
istics  are  indexed  into  a  two  dimensional  array.  Each 
of  the  vehicle  units  is  represented  by  a  row  in  the 
array  and  the  columns  indicate  location  on  the  unit. 

For  calculating  traction,  separate  models  have  been 
developed  for  tracked  and  wheeled  vehicles.  The 
traction  subprogram  (Fig.  17  and  18)  simply  sums  the 
traction  values  for  each  driven  unit  on  the  vehicle. 

Resistance  is  calculated  by  the  motion  resistance 
subprogram  (Fig.  19).  It  is  based  on  eq  31  and  uses 
w  values. 

No  driver  program  has  been  assembled  to  coordi¬ 
nate  the  various  subprograms.  Basically,  the  driver 
program  should  merely  establish  a  user  interface  for 
data  input  and  output  and  call  each  of  the  subprograms 
to  perform  its  function. 


Figure  15.  Inputs  for  wheeled  vehicles. 


Figure  16.  Inputs  for  tracked  vehicles. 
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Figure  1 7.  Traction  subprogram  for  wheeled  vehicles. 


Figure  1 8.  Traction  subprogram  for  tracked,  vehicles. 
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Snow 

type 

Snow 

temp 

'  (*Q 

Air 

temp 

(°C) 

Table  7. 

p 

( y/cm 3) 

Medianicat  properties  of  shallow  snow. 

ca  c  5  <p 

(kPa)  (kPa)  (degrees)  (degrees) 

CO 

(I/cm3) 

New  Snow 

Dry' 

-5  to  -10 

-9 

0.15 

0 

1  to  1.5 

18  to  20 

18  to  20 

0.4  to  0.6 

Dry4 

-10  to -20 

-15 

0.2  to  0.26 

0 

0.5 

22  to  24 

26 

0.7  to  0.8 

Wet5 

-1  to  -3 

-1 

1  to  1.2 

- 

25  to  30 

— 

1  to  1.2 

Dry1 

-3 

-6 

0.13 

0.67  to  1.08 

Dry  b 

-7 

-10 

0.095 

0.972 

Settling  snow 

Powder' 

-3  to  -5 

-8 

0.1  to  0.2 

0 

3.9 

20 

22 

0.75 

Bonded2 

-3  to  -5 

-7 

0.3 

0 

— 

20 

— 

0.6 

Bonded3 

-3  to  -5 

-15 

0.23 

6.8 

25 

25 

0.7 

Powder6 

0 

-2 

0.13 

2.35 

Powderd 

-1 

0 

0.18 

0.68  to  0.72 

Powder* 

-1 

0 

0.18 

1.4.? 

p 

(critical) 

(ti/cm3) 


Old  snow 

Crusted6  -10  -13  0.21 

Crusted6  -5  -10  0.2 

Compacted6  -4  to -10  -1  0.4 

Compacted6  -4  to -10  +6  0.4 

Corn  snow  0  >0  0.6 

'  Slush  0  >0  0.7 

Crusted'  -f  +3  0.36 

Crusted8  -1  +3  0.38 


0 

0.68 

3.0 

0 

0 

0 


Z> 

20 

20 

30 


~  6 

~  4 


1.8S  to  2.11 
2.28 


'Houghton,  30  (an  75. 

2Houghton,  3  Mar  75. 

3Alaska,  Ft.  Greely  and  Ft.  Wainwright,  Feb  72. 
4Houghton,  4  Feb  76, 5  Feb  76. 
sHoughton,  18  Feb  76. 

6Silver  Creek  Bridge,  Feb  73. 

*Alta,  Utah,  27  March  80  (vehicle  test) 
bAlta,  30  March  80  (vehicle  test) 
cAlta,  26  March  80  (vehicle  test) 
dAlta,  1  April  80 
*Alta,  1  April  80  (vehicle  test) 

'Wolf  Creek,  California,  21  March"  80  (vehicle  test) 
8\Volf  Creek,  21  March  80  (vehicle  test) 


Figure  19.  Motion  resistance  subprogram. 
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CONCLUSIONS 

Most  tests  and  measurements  reported  were  con¬ 
ducted  under  conditions  as  they  occurred  during  a 
scheduled  test  period.  This  led  in  most  cases  to  re¬ 
sults  that  were  nongeneric  to  the  terminology  used; 
i.e.  shallow  snow,  packed  snow,  etc.  In  most  cases, 
the  shallow  snow  layers  were  complex  in  nature  and 
varied  considerably  from  site  to  site.  The  scatter  in 
ice  and  packed  snow  data  is  also  indicative  of  varia¬ 
tions  not  accounted  for. 

It  is  therefore  a  reasonable  conclusion  that  there  is 
only  a  small  amount  of  information  available  relative 
to  shallow  surface  layers  in  cold  regions  and  vehicle 
performance. 

Projects  are  now  in  progress  to  establish  a  compre¬ 
hensive  data  bank  over  the  next  three  years  to  correct 
this  deficiency.  Additional  studies  are  planned  to 
establish  mere  accurate  vehicle  performance  models  for 
predicting  mobility  in  winter  environments. 
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A  facsimile. catalog  card  in  library  of  Congress  MARC 
format  is  reproduced  below. 
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